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Preliminary results of one year anisotropy measurement in the energy range 10^^ — 10^* eV as a function of 
energy are presented. The results are compared for two methods of data analysis; the standard one with meteo 
correction approach in use and another one so-called "East minus West" method. Amplitudes and phases of 
anisotropy for three median energies E = 25 TeV, E = 75 TeV and E = 120 TeV are reported. Brief consideration 
of amplitude-phase dependence of anisotropy on energy is expounded. 



1. Introduction 

Experimental results on anisotropy of primary 
cosmic ray flux give important information for 
further theoretical study of the origin and propa- 
gation problems in cosmic ray physics. Together 
with numerical characteristics of anisotropy - am- 
plitude and phase - the energy dependence of 
these characteristic gives additional information 
for better understanding of the above mentioned 
problems. Energy range of primary cosmic rays 
under anisotropy investigation is extremely wide 
- from - 100 GeV up to 10'"^ TeV and higher. At 
1-10 TeV range the measurements are carried 
out with detectors of Cerenkov emission of Exten- 
sive Air Showers or underground muon telescopes 
(Nagoya, Baksan Underground Scintillation Tele- 
scope; MACRO and others). Then, at 100-10000 
TeV range, the EAS-TOP and CASCADE results 
are well known. And finally for 10^^ eV and 
above anisotropy is studied with giant EAS in- 
stallations such as Yakutsk, AGASA and Auger. 
One can find review of anisotropy results in [1] 
and bibliography in [Tl[2l[3]. Some results of 
cosmic ray anisotropy studies with a Small Air 
Shower (SAS) arrays (E ~ 10 TeV) in Northern 
hemisphere were published at seventieth-eightieth 
[HEliniEKH] ■ Amplitude and phase of anisotropy 
were derived from Fourier analysis of counting 
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rate along the Right Ascension coordinate, with- 
out record of SAS arrival direction. The main re- 
sults of those studies were: 1). Compton-Getting 
effect due to revolution of the Earth around the 
Sun in accordance with theoretical prediction [5] 
was observed, 2). Sidereal anisotropy with am- 
plitude ~ 0.06% and phase ^ 1 hour RA was 
measured with high statistical accuracy. When 
there are numerous results for E < 10 TeV and 
E > 100 TeV, the 10-100 TeV range is up to now 
slightly studied. Lack of experimental results at 
this energy range gives no possibility to arrive at 
clear conclusion about dependence of amplitude 
and phase of anisotropy on energy. At 2007, af- 
ter some modernization of Baksan EAS installa- 
tion, we resumed the measurement of anisotropy 
of SAS in the range 10-100 TeV with registration 
of arrival directions on celestial sphere. Below we 
report result of one year registration. 

2. Experimental set-up and Energy Re- 
sponse 

Central part of Baksan SAS array - " Carpet" 

- consists of 400 liquid scintillation detectors, ar- 
ranged in horizontal continuous square geometry 

- Fig.l. Dimension of the detector is 70x70x30 
cm^. There are four outside huts (OH) on square 
diagonals at distance 30 m from the center of Car- 
pet. Each hut contains 18 (3x6)detectors. The 
Carpet is divided into 25 square modules, each 
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module = 16 detectors. Twenty nine (25 mod- 
ules + 4 OH) sum anode pulses (each pulse is 
a sum of anode pulses of individual detectors of 
module) are registered with ADC, TDC and with 
multiple logic unit designed to produce the SAS 
triggers. Individual (each of 29) logic pulse is 
generated if corresponding sum anode pulse ex- 
ceeds level corresponding to energy deposition of 
0.5 relativistic particle. Energy deposition corre- 
sponding to 1 r.p. is equal 50 MeV. Multiplicity of 
individual logic signals together with anode am- 
plitudes are the measure of power of SAS, follow- 
ing the primary particle energy. Arrival direction 
of SAS is calculated using the relative times from 
TDC channels. In the present analysis we have 
used three sorts of triggers: 4-fold coincidence of 
corner modules 4CM, 4-fold coincidence of out- 
side huts 40H, 8-fold coincidence 4CM-f40II. 
Counting rate of coincidences: F(4CM)=1.39 Hz, 
F(4OH)=1.50 Hz, F(4CM-h4OH)=0.72 Hz. An- 
gular resolution: R(4CM)=10°, R(40H)=3.5°, 
R(4CM-h4OH)=3.0°.One year data set was pro- 
cessed, only full days were included in analy- 
sis, useful time after "noise" cuts and elimination 
of "bad" days = 75%. Energy response for three 
sorts of trigger was evaluated by Monte Carlo sim- 
ulations, made on the basis of CORSIKA codes 
(v. 6012, HDPM and Gheisha models). Verti- 
cal proton flux (with differential spectrum 7 = 
- 2.7) uniformly distributed over the circle area 
with r = 60 m was used in simulation. The roof 
of thickness 34 g/cm^above the Carpet was taken 
into account. Evaluated median energies for three 
sorts of trigger are: E (4CM) = 25 TeV, E (40H) 
= 75 TeV, E (4CM-K40H) = 120 TeV. 

3. Processing and analysis of experimental 
information. 

Fourier analysis of resulting day waves in solar, 
sidereal and antisidcrcal time was done to find 
amplitudes and phases of the 1st and 2nd har- 
monics. Resulting day waves were obtained by 
two methods with posterior comparison of meth- 
ods' efficiency: a) traditional approach with mc- 
teo correction, b) differential method " East minus 
West" . The following considerations were used as 
efficiency criteria: 1) Phase of the first solar har- 
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Figure 1. Plan of the Baksan SAS array. 



monic must be close to 6 hours of Local Solar 
Time (Compton-Getting effect), and amplitude 
to be close to 0.034 for our geographic latitude. 
2) Absence of statistically significant antisidereal 
wave. It means the absence of annual modulation 
of solar day wave. Otherwise the false sidereal 
wave inevitably arises as one of two side frequen- 
cies - sidereal and antisidereal. 

3.1. Traditional approach 

Counting rate in cosmic ray variation exper- 
iments depends essentially on variations of me- 
teorological factors - atmosphere pressure and 
temperature. Hence, one need to exclude meteo 
induced variations to get pure initial variations 
of counting rate. Meteorological coefficients for 
three sorts of trigger were derived from multi- 
ple regression analysis. Amplitudes and phases 
of 1st harmonics after meteo correction are pre- 
sented in Table 1. All 2nd harmonic amplitudes 
are inside statistics. It is necessary to note that 
comparatively large statistical error arises from 
short enough data set - only one year of regis- 
tration. And nevertheless, phases of solar waves 
are in good coincidence with expected one from 
Compton-Getting effect - all three are close to 
6 hour of Local Solar Time. At the same time 
extremely large (> 5a) amplitude of 75 TeV an- 
tisidereal wave gives rise to assume that meteo 
correction procedure by itself doesn't guarantee 
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Table 1 



Amplitudes and phases of 1st harmonic after meteo correction. The amplitude errors shown in Tablcl and 
Table2 are statistical ones. 



wave 

1st harmonic 


Solar 


Sidereal 


Antisidereal 


25 TeV amplitude, % 


0.025 ±0.025 


0.065 ±0.025 


0.064 ±0.025 


phase, h 


6.4 ±2.3 


23.5 ± 1.4 


19.1 ± 1.4 


75 TeV amplitude, % 


0.050 ±0.024 


0.088 ±0.024 


0.122 ±0.024 


phase, h 


7.5 ± 1.7 


19.6 ± 1.0 


1.2 ±0.7 


120 TeV amplitude, % 


0.047 ±0.035 


0.059 ±0.035 


0.087 ±0.035 


phase, h 


6.8 ±2.2 


21.1 ±1.7 


22.4 ± 1.5 



total expulsion of extraneous variations. This fact 
arouses doubts about sidereal time result. 

3.2. East minus West method. 

When studying anisotropy we come across not 
only with variations induced by metco factors but 
also with apparatus instabilities. To remove them 
totally or simply to take into account their effect 
is not so easy task. Wit way to solve the problem 
is the 'East minus West" method mentioned for 
the first time in [10] . The method is based on the 
assumption that as meteo factors so apparatus 
instabilities produce equal variations in counting 
rates of showers arriving from East and West di- 
rections. Hence the difference between East-ward 
and West-ward counting rates can eliminate the 
uncontrolled and spurious variations. Sequence of 
differences of counting rates from two directions 
calculated during the day for each fixed time in- 
terval (20-min in our case) results in differential 
day wave. Then integration of resulting (sum of 
all days) differential wave results in primordial 
day wave characterized by amplitude and phase. 
Amplitudes and phases of 1st harmonics resulting 
from East- West method are presented in Table 2. 
Comments to Table 2: 1) Antisidereal waves are 
inside of statistical uncertainties in contrast to 
meteo method. 2) Phases of solar waves are in 
good agreement with Compton-Getting effect. 3) 
There is also a hint that sidereal amplitude de- 
creases with increase of energy - on Fig. 2 points 
of this experiment are connected with dot line. 
Earlier EAS-TOP team also reported compara- 
tively small amplitude [11] for E = 100 TeV. One 



year later Andyrchi result [12] for the same en- 
ergy was presented being in excellent agreement 
with EAS-TOP one. It seems worthwhile to do 
some remarks on possible energy dependence of 
anisotropy. At 30th ICRC the EAS-TOP coUab- 
oration presented [2] results of reprocessed 8-year 
data set using E-W method - point ET07 at Fig. 2. 
We can notice that this result together with pre- 
vious ones for E = 100 TeV really indicate de- 
creasing of amplitude for energy close to 100 TeV. 
And what is more interesting, with further in- 
creasing of energy the amplitude again increases 
up to 0.064 % at E = 400 TeV with drastic change 
in phase [HIS] to 13.6 hour RA from 0.4 hour RA 
for 100 TeV. Apropos, such shift of phase was al- 
ready pronounced in [Tl] for energy 300 TeV 
and higher but with less significance. Phase of 
max counting rate 0.4 hour RA for low energy 
corresponds to direction perpendicular to local 
Orion arm (in out of arm direction) and parallel 
to Galactic plane whereas 13.6 hom- RA points to 
north Galactic Pole inside the local Galactic arm. 
Hence, it is possible to suppose that something 
changes in mechanism producing an anisotropy 
effect in primary cosmic ray flux at energy range 
100 - 300 TeV. Obviously, to arrive at more defi- 
nite conclusion we need acquire more data. 

4. Conclusion 

Posterior analysis demonstrates that "East- 
West" method yields more reliable result on side- 
real anisotropy in comparison with standard me- 
teo correction method. Further accumulation of 
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Tabic 2 

Amplitudes and phases of 1st harmonic after East minus West method. 



wave 

1st harmonic 



Solar 



Sidereal 



Antisidereal 



25 TeV amplitude, % 


0.046 ±0.025 


0.060 ±0.025 


0.045 ±0.025 


phase, h 


6.0 ± 1.9 


21.5± 1.5 


7.9 ± 1.4 


75 TeV amplitude, % 


0.063 ±0.024 


0.043 ±0.024 


0.068 ± 0.024 


phase, h 


6.8 ± 1.4 


22.0 ± 1.9 


5.8 ±0.7 


120 TeV amplitude, % 


0.032 ± 0.035 


0.026 ±0.035 


0.011 ±0.035 


phase, h 


8.3 


0.4 


20.7 
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Figure 2. Amplitude of 1st harmonic of sidereal anisotropy for 10 -100 TeV energy range. 
Ba81 - [Zi], Tib - [13], Ka - [14],Nor - [H], Mus - [5], ET03 - [U], ET07 - [2], An - [12] 



data is necessary to understand better the evolu- 
tion of cosmic ray anisotropy in the energy range 
10 - 100 TeV. Together with resuhs from high 
energy range it will bring some new knowledge 
about origin and propagation of cosmic rays. 
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